Propulsion Mechanisms for Leidenfrost Solids on Ratchet Surfaces 
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We propose a model for the propulsion of Leidenfrost solids on ratchets based on viscous drag 
due to the flow of evaporating vapor. The model assumes pressure-driven flow described by the 
Navier-Stokes equations and is mainly studied in lubrication approximation. A scaling expression 
is derived for the dependence of the propulsive force on geometric parameters of the ratchet surface 
and properties of the sublimating solid. We show that the model results as well as the scaling law 
compare favorably with experiments and are able to reproduce the experimentally observed scaling 
with the size of the solid. 



The propulsion mechanism of self propelling Leiden- 
frost drops on ratchets has been debated since the first 
publication of the phenomenon by Linke et al. in 2006 
PQ . A simpler system that leaves out many of the compli- 
cations due to the deformable nature of the drops [2«5J , 
are Leidenfrost solids. Similar to drops, a platelet of dry 
ice levitating on a cushion of its own vapor over a hot 
ratchet surface shows directed motion [6]. Qualitatively, 
a viscous mechanism was shown to be a likely candidate 
for propulsion in [7] ■ As sublimated vapor is drained side- 
ways along the grooves, perpendicular to the movement 
of the Leidenfrost solid, it is additionally directed down 
the slope of the ratchet towards the deep sections of the 
groove, dragging the levitated solid along. 

In this letter we aim at a more quantitative under- 
standing of the propulsion mechanism for Leidenfrost 
solids. We propose a model based on a continuum de- 
scription for the velocity, pressure and temperature fields 
in the gap between the hot surface and the dry ice. The 
propulsive force is obtained from the viscous drag on the 
surface of the Leidenfrost solid, and we develop a scaling 
expression reflecting how this force depends on geometric 
parameters of the ratchet surface and properties of the 
sublimating solid. This expression is shown to compare 
well with data derived from experiments. 

For our model we assume that the gas velocity it; and 
pressure p as well as the temperature distribution T be- 
tween the hot surface and the platelet of dry ice are gov- 
erned by the steady state incompressible Navier-Stokes 
equations 

pUjdjUi = djlriidiUj + d u % ) - S^p}, (1) 

with diUi = and the energy equation 

pc p u i d i T = d i {Xd i T). (2) 

For simplicity we assume the material parameters to 
be constant. In particular we will use a density of 
p = 1.2kg/m 3 , a viscosity rj = 22uPas, heat capac- 
ity c p — 980J/kg/K and thermal conductivity A = 



29mW/m/K, corresponding to carbon dioxide at a tem- 
perature of 180° C and a pressure of 1 bar. 

Further, we assume the no-slip boundary condition at 
the surfaces of the ratchet and the dry ice. The sub- 
limation rate is determined by equating the heat flux 
normal to the wall J q = —\d n T with the heat sink 
due to sublimation, such that the normal velocity at 
the ice surface becomes u n = Xd n T/(p A_ff subl ), where 
the subscript n designates the normal component and 
AiJ su bi = 598kJ/kg is the latent heat of sublimation. 
We assume the ice surface to be at saturation temper- 
ature, T su bi = — 78.6°C, and the ratchet surface to be 
at a constant temperature of To = 450° C. Despite the 
inhomogeneous sublimation, we also consider that the 
ice surface remains flat. Since the platelet moves across 
the ratchet and hence the surface samples over regions 
of different sublimation rate, this assumption should be 
reasonable. 

The main parameters of interest are the average pres- 
sure p and shear stress r = f zx at the ice surface. The 
first determines which weight of the ice block can be sup- 
ported by the vapor cushion, and the second the viscous 
propulsion force. Using the density 1560 kg/m 3 of dry ice, 
the pressure can be converted into a platelet thickness. 
From the evaporation rate it is apparent that the normal 
velocity scale is governed by Uq = (AAT/A)/( J oA_ff su bi), 
where AT is the temperature difference between the sur- 
faces and A is some length scale of the order of the gap 
width between ratchet and ice, c.f. figure [T] Note in par- 
ticular that the corresponding Reynolds number Re = 
pUoA/r] is independent of this length scale. With the 
above parameters we get Re s» 1. Correspondingly, since 
the forces are assumed to be viscous, the relevant scale 
for the stresses is n ~ r]U a /A ~ ri\AT/(A 2 pAH suh i), 
containing all of the material properties of the vapor and 
the temperature difference. 

To start with, consider only a single groove of width 
L and depth H (c.f. figure [T]). This situation can be 
thought of as a platelet that hovers over many grooves, 
of which one is studied. An ice block of width 2R is 
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FIG. 1. Sketch of the geometry considered. The dry ice 
at temperature Tsubi hovers a distance B above the ratchet 
geometry. The ratchet grooves have periodicity L and depth 
H. In our modeling we consider a single groove of y-extension 
2R with periodicity in x-direction, as well as a disc of radius 
R, covering several grooves. 



FIG. 2. Projected velocity profiles at y = R/2, both xz- 
(top) and yz- velocities (bottom); since u y ^> u x the arrow 
scales are not identical (differing by a factor O(10) at this 
position). R = 5 mm, L — 1 mm, H/L — 0.15, B — 75 pm. 

assumed to levitate at a distance B from the surface. 
Further, we assume the flow, pressure and temperature 
fields to be periodic in the direction normal to the groove 
(x-direction) . At the sides of the ice-block (y = ±R) 
the pressure is ambient and the normal derivatives of 
temperature and velocities vanish. Since R 3> L, H, B, 
the details of this boundary condition do not strongly 
affect the flow inside the groove. 

As an example, we consider the parameters R = 5 mm, 
L = 1 mm, H/L = 0.15 and B — 75 pm. We solve equa- 
tions Q and Q via finite element discretisation as im- 
plemented in Comsol Multiphysics 4.1. The structure of 
the solution is as follows. The main pressure gradient 
is along the groove, underpinning that the groove itself 
acts mainly as a drain for the vapor. The temperature 
profile and the gradients in pressure normal to the y- 
direction become essentially independent of y. Figure [2] 
shows the flow field projected onto the plane normal to 
the groove (arrows), at position y = R/2. Again, this 
distribution is almost independent of y. The flow field 
is essentially as described above: the strong flow along 
the groove is fed by the regions of strong sublimation at 
the narrow gap. The velocity field is almost parabolic, 
albeit slightly squeezed in the xz-projection due to the 
evaporation mass flux at the top wall. 

Since a direct computation of the flow below a realis- 
tically large platelet is beyond our reach, we propose a 
simplified model, assuming - in the spirit of a lubrication 
approximation - a parabolic flow profile everywhere, ne- 



glecting all inertial terms. Defining the average flow rate 
in i direction (i € {x, y}) by integrating over the height 

Mx,y) 

Qi(x,y) = / dzui(x,y,z), (3) 
Jo 

where h(x, y) is the local distance between the ratchet 
and ice surfaces, the Reynolds (lubrication) equation be- 
comes 

dip(x, y) = - h 3^ y \ Qi( x > V)- ( 4 ) 

Mass conservation dictates diQi = u n , with a source term 
due to sublimation equal to the normal velocity at the ice 
surface as defined before. In the spirit of the lubrication 
approximation we assume the local temperature gradient 
to be proportional to l/h(x, y). Taking the divergence of 
the Reynolds equation then yields 

di(h 3 d iP ) = -12(n A 2 )//i. (5) 

We solve equation ^ in the case of a single groove, 
assuming periodicity (p(x,y) = p{x + L,y)), as well as 
in the case of a circular platelet. As boundary condition 
at the edge of the domain f2, i.e. the circumference of 
the platelet, we assume p\aa = and neglect any outlet 
effects beyond that. For the solution a finite element 
discretisation as implemented in Comsol Multiphysics is 
used. In our approximation, the wall shear stress is given 
by t w = - \hVp. 

The lubrication model and the 3D model for a single 
groove give similar shear stress distributions at the upper 
wall. In both cases a net shear stress is generated driving 
the platelet in positive x-direction. To verify the model 
predictions, experiments with levitated dry ice platelets 
have been performed. A quantity that is directly accessi- 
ble experimentally is the acceleration and hence the net 
force on a dry ice platelet of known mass. The stress 
simply is this force divided by ttR 2 . For this, cylinders of 
dry ice with a diameter of 2R = 14 mm ± 0.6mm and 
of different heights are placed on a ratchet heated to 
450 °C. Subsequently, the propelling force is measured 
using the method described in .7 (c.f. supporting in- 
formation) . The geometry parameters of the ratchet are 
L = 1.5 mm and H = 0.25 mm. In figure [3] the exper- 
imentally obtained average wall shear stress is plotted 
vs. the platelet thickness for values between 1 mm and 
10 mm. In the same diagram the corresponding values 
obtained with the 3D model for a single groove and with 
the lubrication model both for a single groove and a circu- 
lar disc are displayed. The model results are unanimously 
compatible with a scaling f zx ~ if-^J of the shear stress 
with the thickness of the platelet, agreeing well with the 
experimental data. However, the different models pre- 
dict slightly different values for the shear stress at a given 
thickness. For the single groove, the lubrication model 
predicts a slightly larger shear stress compared to the 
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FIG. 3. Average shear stress r vs. thickness of ice. Ex- 
perimental data (blue circles), 2D lubrication model (disc: 
dark blue diamonds; single groove: red triangles), 3D simu- 
lation (single groove: green squares). The dotted line shows 
r ~ H? C J. L = 1.5 mm, H = 0.25 mm, R = 7mm. 



more complete 3D model. This is partly due to the ne- 
glect of inertial effects and partly due to the simplified 
temperature profile used in the lubrication model. Com- 
paring the lubrication results obtained for a single groove 
and the more adequate circular disc, the latter predicts 
a larger shear stress. This is due to the fact that for the 
disc a portion of the vapor can also escape to the front 
and back. Therefore the disc hovers closer to the surface, 
resulting in a larger sublimation rate and shear stress r. 
This same effect would also play a role in a more complete 
3D model of a full disc, raising the model results closer to 
the experimental data. We conclude that a model based 
on viscous forces due to the sublimation mass flux not 
only gives the correct order of magnitude for the forces 
on a Leidenfrost solid but is able to reproduce the scaling 
of the forces with respect to the height of the platelet. 

Based on that, we extend the discussion to the scaling 
with respect to the geometric parameters of the ratchet. 
In doing so we assume a constant thickness of the ice 
platelet such that it is supported by an average gas pres- 
sure of p = 100 Pa, which corresponds to a platelet 
of thickness of Hi cc ~ 6.4 mm. The Reynolds equa- 
tion is solved for a circular platelet, and B is varied 
till the average pressure is p. Figure [4] gives an exam- 
ple for the obtained scaling with tana = H/L. Over- 
all, the scaling with geometric parameters is roughly 
t(L,H,R) ~ L 1 /? 3 /? -3 . A similar scaling is obtained 
when considering a single groove in lubrication approxi- 
mation or in the 3D model, albeit with slightly different 
exponents, varying by about 10%. Note that for large 
angles the curve for r vs. tana = H/L deviates from the 
scaling function, and the force increases more rapidly. 
In this region B becomes very small, and eventually the 
pressure below the disc is unable to support its weight, 




H/L = tan(a) 

FIG. 4. Average shear stress r for a disc in lubrication 
approximation as function of the groove angle, tana = H/L 
is varied at fixed L = 1.5 mm and R — 5 mm. B is chosen such 
that the average pressure below the disc becomes p — 100 Pa. 
The dotted line shows r ~ tan 3 a. 
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FIG. 5. Average shear stress r as a function of the platelet 
radius R. A platelet of height H^ ce = 5 mm and radius R = 
7 mm is considered. The red dots are results obtained in the 
2D lubrication model for a disc of corresponding dimensions. 
The dotted line shows r ~ R~ 3 . L — 1.5 mm, H — 0.25 mm. 

i.e. touchdown occurs. In our experiments no movement 
of the discs was observed when increasing the angle to 
tana ~ 0.29, indicating touchdown. Note that a simi- 
lar deviation from the scaling law would be observed for 
small R for the same reason as for large tana. 

To corroborate the findings of the scaling analysis, we 
show in figure [5] the mean stress r measured on a platelet 
as a function of its radius, for a height H lce = 5 mm. Cor- 
responding results obtained in the lubrication approxi- 
mation for a disc are shown in the same figure. The 
stress decreases as the radius increases, and the dashed 
line of slope -3 shows a fair agreement with experiments 
and model results. 

It is interesting to understand this behavior (and more 
generally the scalings found in this study) based on a sim- 
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pie intuitive picture. We divide the flow in a single groove 
into two parts. Firstly, the vapor flows down the incline 
(along the x-axis in fig.l). In the lubrication approxi- 
mation, the velocity of this flow is U\ ~ (A 2 /ry)(APx/ 'L) 
where A is a typical vertical length scale, and APi is de- 
duced from mass flux conservation (c.f. equation (J5J): 
APi ~ r)\ATL 2 /(A A pAH BuU ). The AP X - L 2 scaling 
stems from the quadratic pressure increase with length 
in a channel with constant inflow across permeable walls. 
Secondly, the flow escapes laterally (along the y-axis in 
fig{l]) with a velocity U 2 ~ {H 2 /r])(AP 2 / R). The pres- 
sure associated with this flow balances the weight of the 
platelet (AP2 ~ Pi C c9H lce where p- lce is the dry ice den- 
sity). 

The shear stress arises from the "first" flow, so that it 
scales as rjUi/A. We deduce the length A from mass con- 
servation between the two flows: U\KR ~ U2AL. Even- 
tually, we get an expression for r which is: 



the finite mean free path of gas molecules only play a 
minor role. 



/ pAH suU g 3 p 



1/2 3/2 



H-'^LH 3 

ice 



r]ATX J R 3 ' ^ 

We find that the stress should vary as H^, R~ 3 , H 3 
and L, scalings very close to the ones derived from the 
lubrication approximation (only H^ r J deviates slightly). 
Note that this scaling of equation (pfi is also obtained in 
the lubrication approximation for a single groove in the 
limit of small angles (c.f. supporting information). 

The experimental results and the model both reveal 
a decrease of the propelling force with the radius of the 
platelet (F oc i? _1 ). This behavior stands in contrast 
to what can be observed for Leidenfrost droplets self- 
propelling on ratchets, for which the force was shown 
to increase rapidly with the radius [B]. The systems are 
indeed quite different: a drop on a hot ratchet follows the 
asperities of the texture, so that the characteristic length 
scales involved in the scaling laws above get modified 
- changing the geometry in these systems has a deep 
implication on the resulting force. 

Based on the observations above, we conclude that vis- 
cous drag from pressure driven flow due to sublimation 
seems the most important driving force for self-propelling 
Leidenfrost solids on ratchets. However, despite the ex- 
cellent agreement with scaling predictions our model re- 
sults consistently lie below the measured values, indicat- 
ing that that some detail may not be fully captured yet. 
Also note that these findings do not necessarily carry 
directly over to Leidenfrost drops for which additional 
effects such as deformability or interfacial stresses come 
into play. 

Contrasting this viscous mechanism based on pressure- 
driven flow, it was recently suggested that rarefaction ef- 
fects, in particular thermal creep flow, could play a dom- 
inating role for the propulsion on Leidenfrost ratchets 
[5] . We have investigated these effects using Monte Carlo 
simulations in [9 . For typical geometries used in exper- 
iments it is shown that such rarefaction effects based on 
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SUPPLEMENTARY MATERIAL 



lowest order contribution to the average shear stress as 



Lubrication approximation: single groove 

For a single groove an approximate analytical solution 
of the Reynolds (lubrication) pressure equation ([5| , 

di(h 3 d iP ) = -i2(n A 2 )//i, 

can be obtained by superposing two pressure fields 
p(x,y) — Px{x) +p y (y). This approximation works well 
since, apart from close to the outlet, the x- and y- com- 
ponents of the pressure gradient are almost independent 
of each other. The pressure equation then becomes sep- 
arable and we can obtain a closed form solution. In par- 
ticular, we obtain for the average pressure and the av- 
erage shear stress (using A = L as our length-scale and 
n = ?/AAT/(L 2 pAi7 su bi) as pressure scale): 



p _ AR 2 log Hi /B 
12n £ 2 ~~ ~3~ (Hf - £ 4 ) ' 



1 



12U L 3 



AHBHi 
1 



(7) 



(8) 



H 



2H 2 (Hi+B) G(Hf - B 4 ) 



log Hi/ B, 



where Hi = H + B. Note that since IIq ~ L~ 2 , p be- 
comes independent of L, while f is linear in L in this 
approximation. 

For H <C B equations (|7j and (|8| can be expanded 
to lowest order in H. Solving the expansion of @ for 
B and inserting it into the expansion of yields the 



f = (12IT L 3 ) 



7V3 ( H^p/(12U L 2 ) 



120 



R 



1 Plce9 3 | tj3/2 t „3n-3 



(9) 



where we have used p = pi ce gHi ce . This exactly corre- 
sponds to the scaling found in equation |6| in the main 
text. 



Experimental procedure 

The experiments were performed on the basis of the 
procedure described in [7]. In order to measure the pro- 
pelling force of the ratchet, we push a levitating platelet 
of dry ice of radius R and mass M with some initial ve- 
locity in the direction of negative x (fig. 1). The platelet 
slows down, stops and accelerates in the direction of the 
propulsion. We record its motion with a fast camera (200 
frames per second). We derive the acceleration from the 
trajectory x(t). A constant acceleration x" is measured 
around the point of arrest, from which we deduce the 
propelling force F = Mi". The stress r is defined as 
F/ttR 2 . 

The ratchet is made of brass. Its teeth are machined 
with a milling machine whose head was tilted by 10°. 
The length between two teeth is L = 1.5 mm. We use to 
two different techniques to create the platelet. Either we 
sublimate a small pellet of dry ice of radius R = 7 mm 
to the required height, or we force a block of dry ice into 
holes of different radii drilled into a hot plate. 



